Version: July 13, 2013 

Preprint typeset using DTJjjX style emulateapj v. 11/10/09 


IS THERE A METALLICITY CEILING TO FORM CARBON STARS? - A NOVEL TECHNIQUE REVEALS A 

SCARCITY OF C-STARS IN THE INNER M31 DISK* 

M. L. Boyer 1,2,1 / L. Girardi 1 2 3 , P. Marigo 4 B. F. Williams 5 , B. Aringer 6 , W. Nowotny 6 , P. Rosenfield 5 , C. E. 

Dorman 7 , P. Guhathakurta 7 , J. J. Dalcanton 5 , J. L. Melbourne 8 , K. A. G. Olsen 9 , D. R. Weisz 5 , 

Version: July 13, 2013 

ABSTRACT 

We use medium-band near-infrared (NIR) Hubble Space Telescope WFC3 photometry with model 
NIR spectra of Asymptotic Giant Branch (AGB) stars to develop a new tool for efficiently distinguish- 
ing carbon-rich (C-type) AGB stars from oxygen-rich (M-type) AGB stars in galaxies at the edge of 
and outside the Local Group. We present the results of a test of this method on a region of the inner 
disk of M31, where we find a surprising lack of C stars, contrary to the findings of previous C star 
searches in other regions of M31. We find only 1 candidate C star (plus up to 6 additional, less certain 
C stars candidates), resulting in an extremely low ratio of C to M stars (C/M = (3.3^0^) x 10 -4 ) that 
is 1-2 orders of magnitude lower than other C/M estimates in M31. The low C/M ratio is likely due to 
the high metallicity in this region which impedes stars from achieving C/O > 1 in their atmospheres. 

These observations provide stringent constraints to evolutionary models of metal-rich AGB stars and 
suggest that there is a metallicity threshold above which M stars are unable to make the transition to 
C stars, dramatically affecting AGB mass loss and dust production and, consequently, the observed 
global properties of metal-rich galaxies. 

Subject headings: stars: carbon — stars: AGB and post-AGB — stars: late-type — galaxies: individ- 
ual (M31) 


1. INTRODUCTION 

Asymptotic Giant Branch (AGB) stars play a signifi- 
cant role in galaxies’ observed properties and their evo- 
lution. They are responsible for a major share of galaxy 
luminosity (e.g., Maraston et al. 2006; Marigo et al. 2010; 
Boyer et al. 2011; Melbourne et al. 2012; Melbourne & 
Boyer 2013) and contribute considerably to the chemi- 
cal enrichment of the interstellar medium (e.g., Marigo 
2001; Ventura et al. 2001; Karakas & Lattanzio 2007). 
Despite its widespread importance, the AGB phase re- 
mains among the most uncertain phases of stellar evo- 
lution modeling, leading to the largest uncertainties in 
galaxy stellar population synthesis (Conroy et al. 2009). 

AGB atmospheres have complicated elemental abun- 
dances due primarily to the third dredge-up (3DU) pro- 
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Fig. 1. — The C/M ratio, adapted from Groenewegen (2006a) 
and Battinelli Sz Demers (2005) for nearby galaxies. Data for five 
regions in M31 from Brewer et al. (1995) are plotted with aster- 
isks, with metallicities from Zurita Sz Bresolin (2012). These data 
were derived using optical photometry, so the number of C stars 
is underestimated by ^20-60% (e.g., Boyer et al. 2009). The error 
bars are derived solely from Poisson statistics and do not include 
information on detection biases or photometric uncertainties. The 
solid line is fit through all points, and the dashed line is fit only 
through the Brewer et al. (1995) points. The Davidge et al. (2005) 
point (gray square) and the point from this work (star) are de- 
rived from NIR photometry, with the metallicity range from Saglia 
et al. (2010). We note that the metallicities for M31 are not mea- 
sured directly for the fields included here, but are instead derived 
from estimates of the metallicity gradient across the disk (Zurita 
Sz Bresolin 2012) or the bulge (Saglia et al. 2010). 

cess, wherein they pull newly-synthesized carbon and 
other elements to the surface. The carbon bonds with 
free oxygen to make CO, and the excess oxygen or car- 
bon dictates the molecular and dust chemistry of the 
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Fig. 2. — 3-color image of M31 from Spitzer (Barmby et al. 2006). 24 pm is plotted in red, 5.8 pm in green, and 3.6 pm in blue. The 
PHAT coverage is shown in cyan (Dalcanton et al. 2012b). The coverage discussed here is shown in black and is part of the PHAT survey’s 
“Brick 03”. Four of the regions from Brewer et al. (1995) plotted in Figure 1 are shown in yellow. Region #5 is off the image to the 
southwest. 


star (Marigo & Aringer 2009; Ferrarotti & Gail 2006), 
having drastic consequences for the mass-loss and dust- 
production processes. The ratio (C/M) of C stars (those 
AGB stars with excess free carbon in their atmospheres) 
to M stars (those with excess free oxygen) generally de- 
creases with increasing metallicity (e.g., Rowe et al. 2005; 
Groenewegen 2006a, 2007; Battinelli & Demers 2005). 
Metal-poor stars are more likely to become C-rich be- 
cause less oxygen is available to bind newly dredged-up 
carbon into CO and because the depth of 3DU events 
increases at low metallicity (Karakas et al. 2002). Both 
factors favor the formation of carbon stars after fewer 
and fewer dredge-up events at decreasing metallicity. 
Although this effect is well established on theoretical 
grounds, quantitative predictions of the C/M ratio from 
first principles are not straightforward, given the critical 
dependence of C/M on the star formation history of the 
parent galaxy and the uncertain details of 3DU model- 
ing. Empirically, the relationship between the C/M ratio 
and metallicity is poorly constrained owing to small num- 
ber statistics and other sample biases, especially at high 
metallicity (Fig. 1). 

1.1. Carbon Star Surveys 

Carbon star surveys in nearby dwarf galaxies have 
helped inform evolutionary models at low metallic- 
ity (e.g., Battinelli et al. 2007; Battinelli & Demers 
2009; Groenewegen 2006a, 2007, and references therein; 
Fig. 1). However, the majority of these surveys were con- 
ducted entirely at optical wavelengths (<9000 A) where 
dusty AGB stars are often undetected owing to circum- 
stellar extinction (e.g., Nowotny et al. 2013). In the 
Magellanic Clouds, ~30% of C stars are fainter than 
M\ = —2.5 mag (Boyer et al. 2011), which is the lim- 
iting magnitude of the surveys mentioned above. More 
complete AGB samples can be created using observations 
in the near-infrared (NIR) where circumstellar dust ex- 
tinction is less severe. Several studies have searched for 
carbon stars using NIR photometry in metal-poor dwarf 
galaxies (e.g., Cioni et al. 2006; Sibbons et al. 2012). At 
higher metallicity, NIR searches are much more difficult 
owing either to the distance of the galaxy (in the case 
of M31) or to sample biases and uncertain distance mea- 
surements (in the case of the Milky Way). 

There are two examples of photometrically-complete 
NIR surveys of AGB stars in the inner region of M31: 
Stephens et al. (2003), who used the Hubble Space Tele- 
scope NICMOS camera, and Davidge et al. (2005), who 


used ground-based adaptive optics. Both of these sur- 
veys find examples of stars with very red J — K colors. 
Stephens et al. (2003) identify these stars as long period 
variables, but do not comment on whether they are C or 
M stars. Davidge et al. (2005) identify them as candidate 
C stars based on their NIR colors, suggesting the pres- 
ence of a large C star population. This is in contrast to C 
star surveys in the Milky Way Bulge, which show an al- 
most total lack of intrinsic C stars (e.g., Azzopardi et al. 
1991; Whitelock 1993; Feast 2007; Miszalski et al. 2013). 
Whether or not these red stars in the M31 bulge and 
inner disk are C stars remains unexplored, since narrow- 
band photometry and spectroscopy are difficult both due 
to limited sensitivity and crowding. 

The WFC3 NIR imaging camera on board the Hubble 
Space Telescope has allowed for high-resolution, sensitive 
imaging up to 1.6 /xm, enabling observations of AGB 
stars in galaxies farther than ~5 Mpc for testing and im- 
proving models of AGB mass loss (e.g., Melbourne et al. 
2012). However, it is impossible to distinguish C- and 
M-type AGB stars with the available WFC3 broad-band 
filters alone (Dalcanton et al. 2012a, b). The medium- 
width filters available for WFC3/IR, on the other hand, 
do provide the opportunity to separate AGB subtypes 
(Sect. 2.1). NIR synthetic spectra of AGB stars (Aringer 
et al. 2009, and in preparation) suggest that these filters 
can sample individual molecular features in C-rich AGB 
stars. Therefore, these medium-band filters should pro- 
vide an efficient way of discriminating between C- and 
O-rich AGB stars. 

To test this approach, we targeted a held from Davidge 
et al. (2005) with the WFC3/IR medium-band filters so 
that any candidate C stars detected could be compared 
to those identified via their JHK colors. The Davidge 
et al. (2005) results suggest that a WFC3/IR single held 
should encompass a population of several hundred car- 
bon stars, based on the fraction of AGB stars with red 
J—K and H—K colors in their small 22" held. However, 
our observations reveal the near-absence of C stars. 

We have organized this paper as follows. In Section 2, 
we detail our observations and analysis. Section 3 de- 
scribes our early results regarding the presence of C and 
M-type AGB stars. Finally, the paper is summarized and 
the implications discussed in Section 4. 

2. THE DATA 

Data were collected in a Cycle 20 HST program (GO- 
12862, P.I.: Boyer) and consist of 600 sec total integra- 
tions (300 sec x 2 frames) in each one of the medium- 
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Fig. 3. — Color-magnitude diagrams, plotted as Hess diagrams to show the point-source density. F160W and F110W magnitudes come 
from the PH AT survey (Dalcanton et al. 2012b). Sources with H — K > 0.4 mag from Davidge et al. (2005) are plotted as large green 
circles. The horizontal dotted line is the F160W TRGB from Dalcanton et al. (2012b). The one carbon star candidate identified in Figure 5 
is shown as a red asterisk, and the 6 sources that fall on the edge of the region covered by the carbon star models in Figure 5 are plotted as 
orange diamonds. An additional star on the edge of M star model coverage in Figure 5 is plotted as a blue plus symbol. This star has colors 
consistent with a foreground M dwarf (dashed line, also see Fig. 21 from Dalcanton et al. 2012b), as does one of the C star candidates. 


band WFC3/IR filters: F098M, F127M, F139M, and 
F153M. The 2'3x2'.l field of view is centered at a = 
00 h 43 m 21®7, 6 = +41°21 , 46 , /2, approximately 2 kpc from 
the center of M31 (Fig. 2). Photometry was produced us- 
ing the software package DOLPHOT 2.0 (Dolphin 2000), 
which uses TinyTim point-spread functions (PSFs; Krist 
1995; Hook et al. 2008) to measure stellar fluxes simul- 
taneously in each image. Unreliable measurements were 
culled from the DOLPHOT output using the crowding 
(( C\ ) > 0.25), sharpness (U|5 a| > 1), and signal-to-noise 
{S/N < 3) parameters, resulting in a total of 42, 940 
sources detected in all four filters. 

We performed false-star tests to determine the photo- 
metric completeness and uncertainty, wherein we added 
«50,000 stars (one at a time) using the appropriate PSFs 
into each image and reran the photometry as described 
above. The fake stars were randomly distributed spa- 
tially and mimic the observed color-magnitude distribu- 
tion. The resulting 50% completeness, or limiting mag- 
nitude, is 21.6 < to 50 % < 22.8 mag, with the bluer filters 
showing better completeness. 

To these data, we match the stellar photometric cata- 
logs 12 in the WFC3/IR F110W and F160W filters from 
the Panchromatic Hubble Andromeda Treasury (PHAT) 
program (Dalcanton et al. 2012b). We made no correc- 
tion for interstellar extinction because it is expected to 
be low in these filters and will have minimal effect on 
the separation of C and M stars in F127M— F139M vs. 
F139M— F153M (Sect. 2.1). Barmby et al. (2000) find 
E(B — V) as high as 0.4 mag (or E{ J — K s ) < 0.1 mag; 
Rieke & Lebofsky 1985) for globular clusters near the 
observed field, which can cause a shift of <0.02 mag in 
these colors. Figure 3 shows the color-magnitude dia- 
grams (CMDs) and includes the C star candidates from 

12 http: / /archive. stsci.edu/prepds/phat/ 


Davidge et al. (2005) with H — K > 0.4 (the photome- 
try was redone by Olsen et al. 2006) and candidate C 
stars identified by the medium-band WFC3/IR colors 
(Sect. 2.1). The thin red giant branch indicates mini- 
mal differential extinction across the field. 

2.1. Identifying Carbon Star Candidates 

The location of the change in slope in the luminosity 
function indicates that the tip of the red giant branch 
(TRGB) is M™Gb = 18.45±0.10 mag, or M F153M = 
—6.0 mag for a distance modulus of 24.45 mag (Dalcan- 
ton et al. 2012b). The vast majority of thermally-pulsing 
(TP-)AGB stars should lie above this limit, with excep- 
tions for a subset of stars experiencing minimum flux 
during a thermal pulse cycle and very heavily enshrouded 
stars. We find 3032 AGB stars that are brighter than the 
TRGB; those that are not classified as C stars are con- 
sidered to be M stars. 

To identify C stars, we employ the color combination 
of F127M— F139M vs. F139M— F153M, which is insensi- 
tive to circumstellar and interstellar extinction and takes 
advantage of molecular features present in AGB spec- 
tra. For example, the F153M filter falls entirely within 
a broad absorption feature of C 2 +CN (Fig. 4), with its 
bandhead at 1.4 //m, which is common in C-rich stars. 
O-rich AGB stars of early M subtypes have no signifi- 
cant absorption feature in the range from ~1 to 1.6 /irn, 
whereas those of late M subtypes develop a wide H 2 O ab- 
sorption feature between 1.3 and 1.55 // m which can be 
sampled with the F139M filter. Figure 5b shows the loca- 
tion of model atmospheres for objects of M and C type 
from Aringer et al. (2009, and in preparation) in color 
space, as compared to our WFC3 data. Surprisingly, 
only 1 star resides firmly in the region where C stars are 
expected, and only 6 additional sources fall on the edge 
of the color-color space covered by the C star models 
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(Fig. 5a), far fewer than what is predicted by the relation- 
ship shown in Figure 1 (Sect. 3). It is unclear whether or 
not the 6 stars near the border drawn in Fig. 5b are true 
carbon stars. However, we verify via visual inspection 
that they are real point-sources and are unaffected by 
imaging artifacts and/or crowding. To push any of these 
marginal candidates outside of the shaded region via cor- 
rection for interstellar extinction, E(J — K s ) > 0.8 mag 
is required. This is higher than the expected extinction 
(Sect. 2); however, we note that the shaded region in 
Figure 5 is approximate. 

To the left of the shaded region in Figure 5 and below 
the bulk of the M star population are an additional 8 
sources (0.06 < F127M— F139M < 0.18 mag and 0.13 < 
F139M— F153M < 0.23). Models from Aringer et al. (in 
prep) indicate that these stars may be M dwarfs. TRILE- 
GAL simulations in Dalcanton et al. (2012b) predict that 
a single WFC3/IR field should contain rs 6 foreground M 
dwarf stars. Of these 8 stars, 5 are covered by the cur- 
rently available PHAT catalog, and 3 of these (including 
the one closest to the C star region, marked by a blue 
plus in Fig. 5a) have F110W— F160W colors consistent 
with M dwarf stars (Fig. 3, also see Fig. 21 from Dalcan- 
ton et al. 2012b). Correction for any possible interstellar 
or circumstellar extinction cannot move these stars into 
the C-star region. We thus conclude that none of these 8 
stars are candidate C stars. We also note that one of the 
marginal candidate C stars within the shaded region of 
Figure 5 shows F110W— F160W= 0.57 mag, a color that 
is consistent with M dwarf stars. 

Very dusty carbon stars can be undetected even at 
NIR wavelengths. Reddening vectors in Figure 5 show 
that our method is very insensitive to circumstellar ex- 
tinction. A dust-enshrouded star with E(J — I \ s ) = 
1 mag - for which any color-based C/M classification 
using ground-based JHK photometry would be unreli- 
able - is marginally affected in the WFC3/IR medium- 
band filters, with color excesses FJ(F127M — F139M) and 
£’(F139M — F153M) smaller than 0.2 mag. This has an 
effect only at the edge of the C star region, with the pos- 
sibility that the 6 stars near the C star border are instead 
very dusty M stars. 

In Figure 5, we show only stars above the TRGB. How- 
ever, inclusion of stars up to 1 mag fainter than the 
F153M TRGB reveals no additional C star candidates 
in the shaded region. Therefore, according to the ra- 
diative transfer models of Groenewegen (2006b) and as- 
suming a dust mixture of 70% amorphous carbon + 15% 
SiC, C stars that are undetected here due to circum- 
stellar extinction (those fainter than 1 mag below the 
F153M TRGB) must have mass-loss rates in excess of 
~10 -5 M 0 /yr. In the Magellanic Clouds, <1.5% of C 
stars have mass-loss rates higher than this (assuming so- 
lar gas-to-dust ratio, ip® = 500; Boyer et al. 2012; Riebel 
et al. 2012). Thus, with 1 7 detected C star candidates, 
we expect <0.1 additional dusty C stars in this field. In- 
cluding this in the final C/M ratio results in C/M =(1- 
7.1)/(3025-3032)=(3.3lo°i) x 10 -4 ; with the lower limit 

derived from the uncertainty in the number of M stars. 

3. INTERPRETATION/DISCUSSION 

A few searches for carbon stars in M31 using optical 
photometry with broad and narrow filters have success- 
fully identified AGB stars. Brewer et al. (1995, 1996) car- 
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Fig. 4. — IRTF spectra of Galactic stars (Rayner et al. 2009), 
plotted in gray, for (A) M stars in order of increasing water ab- 
sorption and (B) C stars. Solar metallicity model spectra from 
Aringer et al. (2009, and in preparation) are also plotted for (A) M 
stars (blue/cyan) and (B) C stars (red/orange) with T e ff = 2600 K 
and 3300 K. Shaded regions mark wavelengths affected by tel- 
luric absorption. The model spectra predict all major features 
present in the IRTF spectra. Also shown in black are telluric- 
free spectra of 2 stars in NGC 1850 across the spectral region 
covered by the WFC3/IR G141 grism: panel (A) shows the M 
star 2MASS J05084871-6846220 and panel (B) shows the C star 
2MASS J05085008-6845188. These two stars present 2MASS col- 
ors and magnitudes typical of most C-rich and O-rich AGB stars 
in the LMC and do not show any surprising features in the F139M 
band. 


ried out a photometric narrow-band survey along with a 
spectroscopic follow up of five fields southwest of the cen- 
ter of M31 along the major axis (Fig. 2) and found a to- 
tal of 243 C stars. More recently, Nowotny et al. (2001), 
Battinelli et al. (2003), and Battinelli & Demers (2005) 
conducted additional searches along the same southwest 
axis, using similar techniques. The resulting relationship 
between metallicity and the C/M ratio in M31 is shown 
in Figure 1. 

We can use this relationship to predict the number of 
C stars expected in this single field (not including the 
addition of very dusty, optically obscured C stars that 
are missing from Fig. 1). We assume the surveyed field 
follows the metallicity gradient across the bulge, derived 
by Saglia et al. (2010), who measure the abundances via 
optical long-slit data along 6 position angles, including 
one that passes very near this held (0.05 < [M/H] < 0.2 
at i?M 3 i = 2 kpc). A line fit through all data points in 
Figure 1 thus predicts 40-64 C stars. The range reflects 
both the 3 <r uncertainty in the line fit and the range 
in metallicity. Using only the points from Brewer et al. 
(1995) predicts 11-40 C stars, and the lowest C/M ratio 
shown in Figure 1 yields 57±3 C stars. Here, we investi- 
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gate the discrepancy between these predictions and the 
observed number of C star candidates (Sect. 2.1). 

3.1. Are the missing M31 C stars fainter than our 
detection limit? 

Crowding is smooth across the field of view, at approx- 
imately 3.5 stars/arcsec 2 , so it does not affect the com- 
pleteness at different positions. Based on false star tests 
(Sect. 2), the 50% completeness limit is 21.6-22.8 mag, 
depending on the filter. Since the TRGB is >3 mag 
brighter than the 50% completeness limit at all posi- 
tions, photometric incompleteness has not affected our 
C star sample. Indeed, the false star tests indicate that 
the F153M photometry is 99% complete at the TRGB. 

Circumstellar extinction could also cause carbon stars 
to fall below our detection limit. However, as described 
in Section 2.1, the extreme reddening required indicates 
that we have missed <1 C star due to circumstellar dust. 

3.2. Are we confident that C-star models are correct? 

Our C-M classification (Fig. 5) is heavily based on the 
behavior of the C and M spectral models from Aringer 
et al. (2009, and in preparation), which should there- 
fore be verified. The Infrared Telescope Facility (IRTF) 
database includes NIR spectra of «200 Galactic cool gi- 
ant stars (Rayner et al. 2009). Figure 4 illustrates that 
all observed features are predicted by the model spectra 
from Aringer et al. (2009, and in preparation). However, 
because of significant telluric absorption from water and 
methane from 1.35-1.42 /t m (<20% transmission), spec- 
tral features that cause C and M stars to separate in 
Figure 5 are not observable from the ground. 

As an alternative, we can turn to space-based NIR 
spectral observations. Specifically, AGB stars in the 
LMC cluster NGC 1850 were observed by the HST pro- 
gram 11913 (PI: J. MacKenty) with the WFC3/IR grism 
Gl41, thus providing spectra that avoid the problem 
of telluric absorption. Figure 4b shows the spectra ex- 
tracted for the C star 2MASS J05085008-6845188 ( J - 
K s = 1.65), which reveals no unexpected features in the 
F139M filter that might account for the apparent lack of 
C stars in Figure 5. 

Bright M stars in the same LMC dataset present 
J — K s ~ 1.0 and generally flat spectra across G141, 
as exemplified by the spectrum of 2MASS J05084871- 
6846220 in Figure 4a. These flat spectra are identified 
with M-giant models of T e g > 3300 K, or equivalently 
with early-M subtypes, in which the 1.3-1.55 /mi water 
absorption feature is weak. 

3.3. C stars or late M-giants? 

Our observations show that stars with the reddest 
JHK colors from Davidge et al. (2005) are likely late-M 
giants (i.e., very cool AGB stars, but with O-dominant 
atmospheres) - their colors in Figure 5 point to the 
presence of deep water absorption features, which de- 
press the F139M flux and cause red F139M— F153M and 
blue F127M— F139M colors. The M giant sequence is 
well described by Aringer et al. (in prep) models us- 
ing the Jprgensen et al. (2001) water line list. How- 
ever, difficulty in modeling dynamical atmospheres (pul- 
sation and mass loss) truncates the model sequence at 
F127M— F139M « -0.5 mag. 



F127M-F139M 

Fig. 5. — Color-color diagram comparing the WFC3/IR colors to 
the model predictions for C and M giants. Panel (A) shows the 
observed colors of AGB stars, with red stars from Davidge et al. 
(2005) plotted as green circles and candidate C stars marked as in 
Figure 3. Only stars brighter than the TRGB (Mp 153 M < — 6 mag) 
are included to eliminate most red giant branch stars. Circumstel- 
lar reddening vectors for dust-enshrouded stars are plotted with 
blue (M star) and red (C star) arrows for E(J—K S ) = 1.0 mag (cor- 
responding to a mass-loss rate of «10 — 6 Mq/ yr). These vectors 
are derived from Groenewegen (2006b) radiative transfer models 
applied to cool M and C stellar spectra, with 60% silicate + 40% 
AlOx, and 70% amorphous carbon -j- 15% SiC mixtures, respec- 
tively (see Marigo et al. 2008). In panel (B), AGB star models are 
plotted over the data. The M star models (Aringer et al., in prepa- 
ration) are plotted as blue boxes (using the water line list from 
Jprgensen et al. 2001), with darker colors representing lower surface 
gravity, in the range 2 > logg [cm/s 2 ] > —0.5. The T e ff decreases 
from bottom to top along this sequence, with values in the range 
5000 < T e ff < 2600 K: the kink with the reddest F127— F139M col- 
ors is at T e ff ~ 3300 K. For reference, black pluses show an example 
of the same M star models (logg = 0.5 [cm/s 2 ]) using a different 
water line list (Barber et al. 2006). The C star models (Aringer 
et al. 2009) are plotted as red boxes, with darker colors represent- 
ing higher atmospheric C/O ratios, over the range 1.05 < C/O < 2. 
They span 4000 < T e ff < 2400 K and 0 > log g [cm/s 2 ] > —1, 
with the coolest C star models now being located at the bottom of 
the drawn sequences. All models assume a solar abundance for the 
heavy metals. The shaded region marks the area where C stars are 
expected. 

These stars were identified as C star candidates by 
Davidge et al. (2005) since their red colors ( J — K > 
1.3 mag) are usually markers for C-rich stars in the Mag- 
ellanic Clouds and other nearby galaxies (e.g., Groenewe- 
gen et al. 2009; Nikolaev & Weinberg 2000). However, 
our observations clearly indicate that this is not the case 
in M31. In the metal-poor LMC, the O-rich AGB stars 
are relatively warm and rarely reach spectral subtypes 
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Fig. 6.— TP-AGB Hayashi lines for a model with M = 1.4 Mq at 
varying initial metallicity ([M/H]) and surface C/O ratio, as indi- 
cated. At a given bolometric magnitude, the effective temperature 
is obtained from integrations of complete atmosphere and envelope 
models (Marigo et al. 2013) that include on-the-fly computation of 
both molecular chemistry and opacity (Marigo & Aringer 2009). 

later than M3 (e.g., Elias et al. 1985; Frogel et al. 1990; 
Massey & Olsen 2003). On the other hand, in the metal 
rich M31, O-rich AGB stars fall along a cooler Hayashi 
line (cf. Marigo et al. 2013), which ensures a large frac- 
tion of them appearing as late-M subtypes. The same 
occurs in our own Galaxy and its Bulge, which are known 
to contain M7-M9 giants in copious numbers (e.g., Frogel 
& Whitford 1987; Rayner et al. 2009). This fact is well 
represented in Figure 6, where we see that comparably 
low effective temperatures of a C-rich model with metal- 
licity [M/H] = —0.3 (typical of the LMC), are reached 
by O-rich models with high metallicity (solar and super- 
solar, reasonably suitable for the inner M31 disk). 

There are several examples of known late-M giant stars 
that show very red J — K colors (e.g., Marshall et al. 
2004; van Loon et al. 2005; Javadi et al. 2013), including 
WOHG64, an OH/IR star with T eS ~ 3000 K. How- 
ever, these red M giants are vastly outnumbered by red 
C stars in galaxies such as the Magellanic Clouds and 
M33, where most O-rich AGB stars are warmer and thus 
have early M-type spectral classifications (or even late- 
K). Inspection of optical spectra from Olsen et al. (2011) 
of AGB stars in the Magellanic Clouds indicate that only 
«3% of M stars show H — K > 0.4 mag. The HST 
data indicate that this fraction increases to ~30% in this 
region of M31 (presumably owing to the higher metal- 
licity), assuming that stars that show WFC3/IR colors 
similar to the Davidge et al. (2005) stars will also show 
similar JHK colors. 

3.4. Where are the C stars in this region of MSI ? 

Age effect ? — In old, metal-rich stellar populations like 
the M31 bulge, carbon stars are not expected since the 
efficiency of the 3DU (A) is not sufficient for stars with 
M < 2 M 0 (Karakas et al. 2002). However, from Dor- 
man et al. (2012), we estimate the bulge contributes 
only 28% of the /-band luminosity in this region. The 
/-band shows minimal contamination from short-lived 
bright sources (Melbourne & Boyer 2013), so this frac- 
tion indicates that the bulk of the stellar mass in our 
observed held originates from the younger stellar popu- 
lation of the disk. Models from Courteau et al. (2011) 
fit to the 3.6 [im luminosity profile also indicate that the 
disk population dominates at 2 kpc. 



Fig. 7. — Optical CMD of the observed region, from the PHAT 
program (Dalcanton et al. 2012b). Padova isochrones are plotted 
with A(V) — 0-3 mag for log (age) = 8, 8.5, and 9 (PARSEC, 
version 1.1; Marigo et al. 2008; Girardi et al. 2010). The CMD 
indicates the presence of main sequence turn-off stars at least as 
young as ~300 Myr (or M < 4 Mg). Detailed analysis of this CMD 
is underway (Williams et al. in prep) to quantify the young stellar 
populations in this region, which is not straightforward given the 
level of crowding in the optical HST images. 

PHAT optical CMDs for this region (Fig. 7) also re- 
veal the presence of main sequence turn-off stars younger 
than ~1 Gyr, indicating the presence of AGB stars with 
masses ~3-5 M 0 . With a stellar population of this age, 
we expect a much larger population of C stars based on 
the C/M relationship implied in Figure 1. 

Metallicity effect ? — Brewer et al. (1996) indicate that 
the C/M ratio rapidly decreases towards the inner M31 
disk (down to values of ^0.019 at Rm 31 = 4 kpc) and 
interpret this decrease as being driven by the metallicity 
gradient in M31. Our data at ~2 kpc puts additional 
constraints in the way the C/M ratio decreases towards 
the M31 center. Within 2 kpc of M31’s center, Saglia 
et al. (2010) find that the metallicity spans 0 < [M/H] < 
0.5, with [M/H] a; 0.1 ± 0.1 near the position of our 
observations (see their Fig. 12). If our very low C/M 
ratio is caused by the high metallicity at this radius, it 
supports the presence of a metallicity threshold above 
which C stars simply do not form. 

A reduction in the C /M ratio with metallicity is indeed 
indicated by all modern TP-AGB models (Karakas et al. 
2002; Marigo et al. 2013) as a result of two effects: the 
larger amount of carbon that needs to be dredged-up to 
make the C/0>1 transition, and the fact that the 3DU 
starts later (at higher luminosities) and is less efficient at 
increasing metallicity. Figure 8 shows the C/M ratio ac- 
cording to the new TP-AGB models from Marigo et al. 
(2013). For each combination of age (or initial mass) 
and initial metallicity, the quantity C/M is obtained as 
tq /tm , he. the ratio of the lifetime spent by the star in 
the C-rich mode (C/0> 1) over the one spent in the O- 
rich mode (C/0< 1). To avoid contamination from the 
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Fig. 8. — Map of the C/M ratio in the age-metallicity plane, 
according to the new TP-AGB models computed by Marigo et al. 
(2013). The lack of the C-star phase (C/M = 0) corresponds to 
the darkest blue region, which extends towards younger ages at 
increasing metallicity. Above some critical value of [M/H], models 
predict C/M= 0 over the entire age range. According to Saglia 
et al. (2010), the metallicity of the region observed here is [M/H] « 
0.1 zb 0.1. 

RGB phase, both tq and tm are evaluated for bolometric 
magnitudes Mboi < —3.6 mag. The region of C/M = 0 
(dark blue in Fig. 8) extends to younger ages at increas- 
ing metallicity, encompassing the whole age range for 
[M/H] > 0.4. At these metallicities, more massive TP- 
AGB stars (ages < 1 Gyr) go through a few dredge-up 
events, but the termination of the TP-AGB phase (deter- 
mined by the onset of the super-wind mass loss) occurs 
before C/O can exceed unity. The models from Marigo 
et al. (2013) also indicate that, while the trend depicted 
in Figure 8 has a general validity, the exact upper limit 
in metallicity for the formation of C stars is sensitive to 
model details, mainly related to the 3DU and mass loss. 
In this context, observations like the ones we present in 
this work are of paramount importance to set stringent 
constraints on the evolution of TP-AGB stars in the still 
largely unexplored regime of high metallicity. 

4. CONCLUSIONS 

We have conducted a search for carbon-rich AGB stars 
in the inner disk of M31 using the WFC3/IR medium- 
band filters on board HST. We find that: 
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